RESEARCH ARTICLE 



Livestock Origin for a Human Pandemic Clone of Community- 
Associated Methicillin-Resistant Staphylococcus aureus 

Laura E. Spoor, 3 Paul R. McAdam, 3 Lucy A. Weinert, b Andrew Rambaut, c Henrik Hasman, d Frank M. Aarestrup, d Angela M. Kearns, e 
Anders R. Larsen/ Robert L. Skov, f J. Ross Fitzgerald 3 

The Roslin Institute and Edinburgh Infectious Diseases, University of Edinburgh, Easter Bush, Midlothian, United Kingdom"; University of Cambridge, Department of 
Veterinary Medicine, Cambridge, United Kingdom 13 ; Institute of Evolutionary Biology, Ashworth Laboratories, University of Edinburgh, Edinburgh, United Kingdoms 
National Food institute, Technical University of Denmark (DTU), Lyngby, Denmark d ; Microbiology Services, Colindale, Health Protection Agency, London, United 
Kingdom 0 ; Department of Antimicrobial Surveillance and Research, Statens Serum Institute, Copenhagen, Denmark' 

ABSTRACT The importance of livestock as a source of bacterial pathogens with the potential for epidemic spread in human popu- 
lations is unclear. In recent years, there has been a global increase in community-associated methicillin-resistant Staphylococcus 
aureus (CA-MRSA) infections of healthy humans, but an understanding of the different evolutionary origins of CA-MRSA clones 
and the basis for their recent expansion is lacking. Here, using a high-resolution phylogenetic approach, we report the discovery 
of two emergent clones of human epidemic CA-MRSA which resulted from independent livestock-to-human host jumps by the 
major bovine S. aureus complex, CC97. Of note, one of the new clones was isolated from human infections on four continents, 
demonstrating its global dissemination since the host jump occurred over 40 years ago. The emergence of both human 5. aureus 
clones coincided with the independent acquisition of mobile genetic elements encoding antimicrobial resistance and human- 
specific mediators of immune evasion, consistent with an important role for these genetic events in the capacity to survive and 
transmit among human populations. In conclusion, we provide evidence that livestock represent a reservoir for the emergence of 
new human-pathogenic S. aureus clones with the capacity for pandemic spread. These findings have major public health impli- 
cations highlighting the importance of surveillance for early identification of emergent clones and improved transmission con- 
trol measures at the human-livestock interface. 

IMPORTANCE Animals are the major source of new pathogens affecting humans. However, the potential for pathogenic bacteria 
that originally were found in animals to switch hosts and become widely established in human populations is not clear. Here, we 
report the discovery of emergent clones of methicillin-resistant Staphylococcus aureus (MRSA) that originated in livestock and 
switched to humans, followed by host-adaptive evolution and epidemic spread in global human populations. Our findings dem- 
onstrate that livestock can act as a reservoir for the emergence of new human bacterial clones with potential for pandemic 
spread, highlighting the potential role of surveillance and biosecurity measures in the agricultural setting for preventing the 
emergence of new human pathogens. 
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Emerging infections represent a major challenge to human and 
veterinary medicine. The great majority of new pathogens of 
humans originate in animal populations, but most are associated 
with episodic zoonotic infections that do not have the capacity to 
transmit to other individuals ( 1 ) . Staphylococcus aureus is a major 
pathogen responsible for considerable human morbidity and 
mortality on a global scale (2). S. aureus is also a leading cause of 
infections of livestock such as cows and is a major economic bur- 
den on the global dairy industry (3). The results of population 
genetics studies have shown that most strains of S. aureus are host 
specific, indicating a low frequency of cross-species transmission 
(3). However, recent studies employing multilocus sequence typ- 
ing (MLST) and whole-genome sequencing (WGS) have identi- 
fied several S. aureus sequence types (ST) that are associated with 
multiple host species, implying either zoonotic transmission or a 



recent common ancestor. For example, the poultry-associated se- 
quence type 5 (ST5) and livestock-associated, methicillin- 
resistant S. aureus (LA-MRSA) ST398 clones are descended from 
bacteria that recently made human-to-animal host jumps (4, 5). 
Importantly, LA-MRSA ST398 strains have acquired antibiotic 
resistance due to selective pressures in the pig farming industry 
and have the capacity to cause severe zoonotic infections of hu- 
mans in contact with pigs on pig farms (5). However, carriage of 
LA-MRSA ST398 by pig farmers is transient, and LA-MRSA 
ST398 does not appear to be readily transmissible between hu- 
mans, probably due to the loss of expression of proteins associated 
with the cell wall that are required for human host colonization 
and transmission (6). Loss of function of superfluous genes is a 
common feature of bacteria undergoing niche adaptation (7), 
which is likely to attenuate the capacity to infect the ancestral host 
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species (6). Recently, Harrison et al. demonstrated that ST130 
MRSA isolates with the mecC gene have spread from cows to hu- 
mans resulting in clinical infections (8). However, the potential 
for animal-specialized strains of S. aureus to successfully cross the 
human species barrier and become epidemic in human popula- 
tions is unclear. Previously, we reported that the clonal complex 
59 (CC59) S. aureus clone that is endemic in Taiwan and has 
spread to other parts of the world may have originated in livestock 
about 500 years ago, and this has recently been corroborated by 
Shepheard et al. (9). However, the limited resolution of MLST has 
precluded rigorous examination of the occurrence of more-recent 
livestock- to-human host jump events leading to the emergence of 
new epidemic S. aureus clones (10). 

S. aureus strains belonging to MLST CC97 are a leading cause 
of bovine mastitis in Europe, Asia, and North and South America 
(11, 12). Less commonly, CC97 has also been reported to cause 
infections of small ruminants, pigs, and humans (11, 13, 14). Im- 
portantly, there have been increasing reports of human infections 
caused by CC97 isolates in Europe, North and South America, 
Africa, and Asia (13, 15-18). However, the origin of human CC97 
strains and their relatedness to livestock-associated CC97 strains 
are unknown. Here, we investigate the evolutionary history of the 
CC97 S. aureus lineage and identify clones that are epidemic in 
human populations that evolved through host jumps from cows 
followed by human host adaptation. Our findings highlight cows 
as a potential reservoir for the emergence of new clones with the 
capacity for pandemic spread in humans. 

RESULTS AND DISCUSSION 

CC97 is an emerging cause of human infections. The earliest re- 
port in the literature of a CC97 isolate is the archetypal bovine 
reference strain of S. aureus, Newbould 305 (NCIMB 702892), 
isolated from a case of mastitis in Canada in 1958 (19), whereas the 
earliest report of human CC97 isolates relates to strains isolated 
38 years later in 1996 (20). However, human CC97 S. aureus bac- 
teria have since been identified in at least 35 countries (see Ta- 
ble S2 in the supplemental material). In order to investigate the 
possibility that infections due to CC97 strains are increasing in 
prevalence in human populations, we determined the number of 
CC97-associated strains isolated from MRSA and bacteremia in- 
fections in Denmark between 2007 and 2011. Uniquely, in Den- 
mark, submission of all MRSA to the Danish National MRSA Ref- 
erence Laboratory, Statens Serum Institut, has been mandatory 
since November 2006, and all MRSA and bacteremia S. aureus 
isolates submitted since 2007 have been genotyped by staphylo- 
coccal protein A (spa) typing. In Denmark, cases of MRSA caused 
by CC97-related spa types increased from a total of 2 in 2007 to 22 
in 2011, equivalent to an 11 -fold increase in 5 years. This repre- 
sents a significant increase in prevalence from 0.3% to 1.7% of 
total annual MRSA in Denmark since 2007 (P < 0.01). Al- 
though surveillance systems for other countries are currently 
inadequate to identify trends in S. aureus genotype associated 
with human infections, a recent study of community-associated 
S. aureus isolated from humans in 16 countries in Europe re- 
ported that 3% of community-associated MRSA (CA-MRSA) and 
8% of community-associated methicillin-sensitive S. aureus (CA- 
MSSA) were ST97 isolates (21). Taken together, these data suggest 
that CC97 is an emerging cause of human infections. 

Human ST97 S. aureus clones originated from bovine-to- 
human host jumps. In order to investigate the evolutionary his- 



tory of the CC97 clone, we obtained 220 CC97 S. aureus isolates of 
human, bovine, porcine, and caprine origin, isolated in 18 coun- 
tries on four continents between 1956 and 2012. For whole- 
genome sequencing, we selected 43 CC97 S. aureus isolates (in- 
cluding 16 MRSA) which broadly represented the breadth of host, 
geographic, and temporal variation identified among strains re- 
ported in the literature, in addition to a single isolate of the closely 
related ST28 as a human-specific outgroup. The core genome of 
the 43 CC97 isolates consisted of 2,079,972 bp which included 
5,425 high-quality single nucleotide polymorphisms (SNP). Phy- 
logenetic reconstruction of the CC97 lineage was carried out by 
using the BEAST program (named BEAST for Bayesian evolution- 
ary analysis by sampling trees) (22) (Fig. 1). The consistency index 
for parsimonious sites predicted using PAUP*4.0blO (23) indi- 
cated a very low predicted frequency of homoplasic events (con- 
sistency index [CI] = 0.95), implying an unconfiicted phyloge- 
netic signal. Furthermore, a maximum likelihood tree constructed 
from the same core genome alignment resulted in a largely similar 
phylogeny with strong bootstrap support at most nodes (see 
Fig. SI in the supplemental material). Overall, the high-resolution 
phylogenetic tree resolves S. aureus CC97 into distinct host- 
associated clades (Fig. 1). Of note, there is considerable genetic 
diversity among CC97 isolates of livestock origin, which is indi- 
cated by numerous deep branches in the phylogenetic tree, com- 
pared to human isolates which are restricted to two distinct clades 
of closely related isolates (Fig. 1). The majority of livestock- 
associated CC97 isolates lie basal to the human clades, and a 
continuous-time Markov model with host association as a discrete 
trait (10, 36) provided strong support that the most recent com- 
mon ancestor of each human clade was of bovine origin (see Ta- 
ble S3 in the supplemental material). These data indicate that 
CC97 isolates circulating among human populations are the result 
of livestock- to-human host jumps which have occurred on at least 
2 independent occasions. Human CC97 clade A consists of isolates 
of both MSSA and MRSA from 12 different countries on four 
different continents, indicating its global dissemination (Fig. 1). 
In contrast, human CC97 clade B is represented by a single MSSA 
isolate from a Danish patient in 1980 and several MRSA bacteria 
isolated within the last 10 years in the Midlands region of the 
United Kingdom, indicating a more limited geographic distribu- 
tion among the strains sampled (13). In order to determine the 
time frame of the livestock- to-human host jump events, we deter- 
mined mutation rates for the CC97 lineage, allowing for variation 
in rates associated with different clades. Initially, an uncorrelated 
lognormal relaxed molecular clock model was used to determine a 
mutation rate for the livestock strains only of 1.53 X 10 -6 nucle- 
otide substitutions per site per year (95% highest posterior densi- 
ties [HPDs] of 1.35 X 10" 6 to 1.72 X 10" 6 ). This rate was then 
fixed for livestock strains, and a local rate clock model was applied 
to each of the human clades, resulting in estimates of 9.58 X 10~ 7 
(7.80 X 10" 7 to 1.15 X 10" 6 ) for human clade A and 1.29 X 10~ 6 
(1.05 X 10" 6 to 1.54 X 10" 6 ) for human clade B. These data 
indicate that closely related S. aureus clones can have different 
mutation rates which may reflect variations in selective pressures 
encountered in different ecological niches. The Bayesian analysis 
resulted in estimates of host jump events which occurred between 
1894 and 1977 for human CC97 clade A and between 1938 and 
1966 for human CC97 clade B (Fig. 1; Table S3). In addition, we 
estimated the date of the most-recent human ancestor (MRCA) of 
CC97 and the human outgroup ST28 to be AD 784 (BC 325 to 
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FIG 1 Identification of human epidemic S. aureus clones descended from bacteria that made livestock-to -human host jumps. The Bayesian phylogenetic 
reconstruction of the CC97 lineage is shown. The tree is based on core genome alignment with branches color coded according to the host species association 
(blue, bovine; green, porcine; orange, caprine; red, human) and date and country of origin of each isolate indicated. The presence or absence ( — ) of /3-toxin phage 
IEC variants B and E and SCCmec type IV or V is indicated by the appropriate letter, and the presence of 5. aureus pathogenicity island (SaPI) -encoded vwb or 
phage-encoded lukM/lukF is denoted by black and gray circles, respectively. The branch lengths are scaled according to the time scale bar (years) and the posterior 
probability values are indicated at each node. Clades A and B are shown. 



AD 1460), which lies within the credibility intervals of the previ- 
ous estimate of the most recent minimum date for the host jump 
event using MLST data (Table S3) (10). 

The results of previous studies employing MLST have sug- 
gested that livestock-associated strains of S. aureus originated 
from human ancestral strains through human-to-animal host 
jumps leading to host-adapted clones specialized for livestock 
with occasional host jumps back into humans ( 10). In particular, 
the CC59 S. aureus clone that is endemic in human populations in 
Taiwan may have originated in livestock about 500 years ago (10). 
Here, we provide an example of recently emerged clones of human 
S. aureus that evolved through independent host jumps from live- 
stock. These data demonstrate that livestock represent a potential 
reservoir of pathogenic bacteria that may cross the species barrier 
and spread among global human populations. 

Methicillin resistance was acquired by human CC97 clones 
subsequent to the host jumps from cows. Antibiotic susceptibil- 
ity testing of CC97 isolates revealed that 7 of 17 bovine isolates 
were sensitive to all antimicrobial agents tested, with a further 6 
resistant to only a single agent, demonstrating the low prevalence 
of antimicrobial resistance among isolates of the leading bovine 
clone of S. aureus (see Fig. S2 in the supplemental material). In 
contrast, 20 of 23 human isolates were resistant to at least one 
antimicrobial agent, with some strains resistant to j3-lactam anti- 
microbials, lincosamides, erythromycin, and trimethoprim. Im- 



portantly, methicillin resistance is a characteristic of the human 
CC97 clades, with two distinct staphylococcal cassette chromo- 
some mec element (SCCmec) types, types IV and V, which are 
associated with human CC97 clades A and B (Fig. 1). Of note, 
none of the isolates examined contained the novel mecC allele 
previously found among several bovine MRSA clones responsible 
for episodes of human zoonotic infection (8). In contrast, all of the 
bovine S. aureus isolates examined in the study are methicillin 
sensitive, suggesting that resistance was acquired after the host 
jump from cows to humans, presumably as a result of selective 
pressures imposed by prescription of antibiotics for treating hu- 
man infections. Consistent with this, the earliest human isolate 
identified in cladeB (isolated in Denmark in 1980) was sensitive to 
all antibiotics tested, whereas all other isolates of clade B (isolated 
since 2005) were resistant to multiple antimicrobial classes (see 
Fig. S2 in the supplemental material). Recently, it was demon- 
strated that methicillin and tetracycline resistance was likely ac- 
quired by LA-MRSA ST398 strains by antibiotic selective pres- 
sures encountered within the pig farming industry (5). Of note, 
the two ST97 isolates from pigs in the current study were resistant 
to both tetracycline and methicillin, in addition to ciprofloxacin 
(Fig. S2). Overall, the antimicrobial susceptibility profiles of the 
human and pig CC97 isolates demonstrated resistance to a much 
greater number of antimicrobials than the bovine CC97 S. aureus 
(Fig. S2). These data imply that the dairy industry does not 
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strongly promote the emergence of antibiotic-resistant S. aureus, 
in spite of the widespread use of antibiotics for treating bovine 
mastitis. We speculate that this reflects the ability of bovine strains 
to invade and survive within bovine mammary epithelial cells, a 
niche which may have a markedly reduced antibiotic selective 
pressure (24). 

The content of MGE among CC97 isolates correlates with 
host species. Previous studies have highlighted the importance of 
mobile genetic elements (MGE) in adaptation of S. aureus to dis- 
tinct host species (3). In the current study, of the 20 livestock 
isolates, 3 and 8 isolates contained MGE encoding the LukM/F 
leukotoxin and the von Willebrand binding protein (vWbp), re- 
spectively, both of which have been demonstrated to have rumi- 
nant host-specific activity (3) (Fig. 1). In contrast, none of the 
human isolates contained the livestock- associated MGE. The fam- 
ily of j8-toxin-converting phages (9Sa3) associated with human 
S. aureus typically contain an immune evasion cluster (IEC) of 
genes encoding secreted proteins such as staphylokinase (sak), 
staphylococcal complement inhibitor (sen), and chemotaxis in- 
hibitory protein of S. aureus {dip) which contribute to immune 
evasion in a human host-specific manner (25). Consistent with 
previous reports, in the current study, we found that 19 of 23 
human isolates and none of 19 bovine or pig isolates contained a 
cpSa3 with an IEC, demonstrating a strong correlation with a hu- 
man host association (25). Of the 16 human CC97 clade A isolates, 
14 had an IEC containing genes sak and sen (IEC type E), and of 
the 8 human CC97 clade B isolates, 4 had an IEC containing sak, 
chp, and sen (IEC type B) and one isolate had an IEC type E (Fig. 1 ) 
(25). Of note, the single goat strain that cosegregates with human 
clade A contains a tpSa3 phage which is highly similar to that of the 
human strains in clade A, implying that it is a human contaminant 
or the result of a very recent human-to-goat transmission event. In 
addition, the arginine catabolic mobile element (ACME) is a char- 
acteristic of some CA-MRSA clones, including the highly success- 
ful USA300 clone which may contribute to enhanced survival dur- 
ing community- associated infections (26). Of note, all seven 
United Kingdom isolates of the human CC97 clade B contained 
ACME linked to SCCmec type V (13). Taken together, the distri- 
bution of MGE among CC97 isolates reveals horizontal gene ac- 
quisition events which correlate with adaptation to different host- 
associated ecological niches. 

Concluding comments. The recent global increase in CA- 
MRSA infections reflects the expansion of an array of S. aureus 
clones with distinct evolutionary histories. Here, we demonstrate 
that livestock is one potential reservoir of pathogenic bacteria with 
the capacity to cross the species barrier, undergo host-adaptive 
evolution, and become established in global human populations. 
Furthermore, the data suggest that a limited number of genetic 
events may be sufficient to transform an S. aureus strain which has 
coevolved with bovine hosts over several thousand years into a 
successful human epidemic lineage. 

The importance of hygiene in prevention of hospital transmis- 
sion of nosocomial pathogens such as MRSA is widely appreci- 
ated, and the recent reduction in hospital MRSA infections is 
likely due in part to improved hygiene measures for controlling 
transmission (27). Improved biosecurity and hygiene control 
measures which prevent the spread of bacterial flora between live- 
stock and human hosts may limit opportunities for successful 
livestock-to-human transmission. Furthermore, regular surveil- 
lance of the microbiota in livestock and humans may facilitate the 



early identification of emergent clones with the capacity to trans- 
mit and cause disease among human populations. 

MATERIALS AND METHODS 

Bacterial isolates. A total of 220 S. aureus isolates of CC97 from bovine, 
human, porcine, and caprine hosts isolated between 1956 and 2012 in 18 
different countries on four continents were obtained. For whole-genome 
sequencing, a total of 43 CC97 strains were selected to represent the 
breadth of host, clinical, spatial, and temporal variation (see Table SI in 
the supplemental material). Strains were grown for 16 h on tryptic soy 
agar (TSA) at 37°C or in tryptic soy broth (TSB) at 37°C with shaking at 
200 rpm. Genomic DNA was isolated using the PurElute bacterial 
genomic kit (Edge BioSystems, MD), with an amended protocol as previ- 
ously described (4). MLST was conducted to confirm the sequence type 
(ST) prior to whole-genome sequencing using methods described previ- 
ously (28). 

Genome sequencing, mapping assembly, and SNP calling. Paired- 
end Illumina sequencing of bacterial strains was carried out on an Illu- 
mina genome analyzer IIx or Miseq, following standard Illumina proto- 
cols. Nucleotide distribution and quality scores of raw reads were assessed 
and filtered for quality using the FASTX toolkit (http://hannonlab.cshl 
.edu/fastx_toolkit/index.html) (genome assembly metrics are provided in 
Table S4 in the supplemental material). High-quality reads were aligned 
against the reference genomes of S. aureus strains MW2 (accession num- 
ber NC_003923) orMu50 (accession number NC_002758.2), variant sites 
were detected, and consensus sequences were called as described previ- 
ously (29), with extended base alignment quality (BAQ) computation 
employed in SAMtools v.0.1.16 (30). The core genome was defined as all 
nucleotide sites shared by all isolates, and putative recombinant regions 
were removed after breakpoint detection with the suite of programs in- 
cluded in the Recombination Detection Program v4.13 (RDP) (31). For 
comparisons of gene content, de novo assemblies of short reads was 
carried out using Velvet vl.0.15 (32) and the VelvetOptimiser.pl script 
within VelvetOptimiser v-2.1.7 (http://bioinformatics.net.au/software 
. velvetoptimiser.shtml) . 

Phylogenetic analysis. Core genome and in-frame protein-coding 
sequences were extracted from the genome consensus sequences to 
construct alignments using custom scripts. The maximum likelihood 
phylogeny was reconstructed using RAxML-7.2.6 (33) implementing a 
generalized time reversible (GTR) model with gamma correction for rate 
heterogeneity, and 1,000 bootstrap replicates (see Fig. SI in the supple- 
mental material). Bayesian phylogenetic analysis was conducted using 
BEAST vl.7.1 (34) implementing the Hasegawa-Kishino-Yano model of 
sequence evolution with gamma correction for rate heterogeneity. The 
rates were calculated by the dated tip method using a three-rate local clock 
model, with the livestock rate constrained to 1.53 X 10~ 6 nucleotide 
substitutions per site per year as determined using an uncorrelated relaxed 
molecular clock model with a constant coalescent prior (35). For the date 
of the MRCA with the human outgroup ST28, an in-frame coding se- 
quence (CDS) alignment of CC97 and ST28 sequences was employed 
which allowed the estimation of mutation rate over all sites and third 
codon positions only. In each case, the mutation rates were comparable 
and resulted in similar estimates of the host jump date (Table S3). Ances- 
tral host states were predicted using a model of discrete trait evolution 
(36), adapting the method by substituting geographic location with host 
state as described previously (10). The consistency index for parsimony 
informative sites was determined using PAUP v4.0bl0 (23). 

Accessory genome analysis. De novo contigs and reference assembly 
sequences were interrogated for specific MGE using BLASTn (37), and 
MGE sequences were aligned using Mauve v2.3.1 (38). Immune evasion 
cluster (IEC) type was defined as described previously (25). 

Antimicrobial sensitivity testing. All isolates in the study were tested 
using the Vitek 2 system (AST-P620 card) (bioMerieux United Kingdom 
Limited, Basingstoke, United Kingdom) with a panel of antimicrobial 
agents, including cefoxitin, benzylpenicillin, oxacillin, gentamicin, cipro- 
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floxacin, inducible clindamycin resistance, erythromycin, clindamycin, 
linezolid, daptomycin, teicoplanin, vancomycin, tetracycline, tigecycline, 
nitrofurantoin, fusidic acid, mupirocin, chloramphenicol, rifampicin, 
and trimethoprim. Results for antibiotic susceptibility were interpreted 
using standards defined by the Clinical and Laboratory Standards Insti- 
tute (CLSI) (39). 

Statistical analysis. For analysis of the increase in CC97-associated 
infections in Denmark, the Cochran-Armitage trend test was employed 
(GraphPad Prism 5.0). 

Nucleotide sequence accession numbers. The Illumina sequences 
generated and used in this study are deposited and available in the Se- 
quence Read Archive (SRA) (http://www.ncbi.nlm.nih.gov/sra) under the 
study accession number PRJEB141 1, located at http://www.ebi.ac.uk/ena 
/data/view/PRJEB141 1. The S. aureus isolates are under sample accession 
numbers ERS212248 to ERS212287 and ERS249844 to ERS249847. 

SUPPLEMENTAL MATERIAL 

Supplemental material for this article may be found at http://mbio.asm.org 
/lookup/suppl/doi: 10.11 28/mBio.00356- 1 3/-/DCSupplemental. 

Figure SI, PDF file, 0.1 MB. 

Figure S2, PDF file, 0.2 MB. 

Table SI, DOCX file, 0.1 MB. 

Table S2, DOCX file, 0.1 MB. 

Table S3, DOC file, 0.1MB. 

Table S4, PDF file, 0.2 MB. 
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